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ABSTRACT 


The  propagation  of  an  initially  sharp  plane  pressure  pulse  through 
a  linear  elastic  composite  medium  is  analysed.  Wave  front  and  ray  theory 
analogous  to  geometrical  optics  is  shown  to  determine  the  change  in 
shape  of  the  leading  wave  front  and  also  the  stresses  immediately  be¬ 
hind  it.  For  certain  circumstances  the  stress  amplitudes  on  this  front, 
or  the  corresponding  tensile  stresses  on  its  reflection  at  the  free 
back  surface  of  a  slab,  may  be  critical  in  design.  Examples  are  presen¬ 
ted  of  an  initially  sharp  plane  pressure  pulse  transmitted  through  an 
elastic  circular  cylinder  and  an  elastic  spherical  inclusion.  The 
method  can  be  applied  to  more  general  composite  configurations,  and 
can  be  extended  to  determine  the  stress  gradient  behind  the  front. 


Fbr  the  latter,  general  formulae  are  derived  by  which  the  reflection 
and  transmission  coefficients  can  be  determined  for  the  stress  gradient 


and  the  higher  order  derivatives  at  an  arbitr 
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1. 


Introduction. 


The  damage  which  an  initially  plane  compressive  stress  pulse  can 
create  when  passing  through  a  slab  of  composite  material  will  depend 
on  the  ■ju.wenuat jun  ut‘  the  pulse  due  to  absorption  and  scattering,  and 
this  influence  on  the  tensile  stresses  associated  with  reflection  at 
the  back  free  surface  of  the  slab  which  can  cause  fracture  and  spalling. 
Few  analytical  investigations  of  such  situations  appear  to  have  been 
carried  out,  and  we  therefore  present  a  study  which  provides  information 
on  certain  aspects  of  the  phenomenon.  In  particular,  we  consider  the 
variation  in  shape  and  stress  magnitude  of  the  wave-front  as  it  traverses, 
for  instance,  a  regular  array  of  parallel  cylindrical  rods  embedded 
parallel  to  the  wave  front,  as  depicted  in  Fig.  la.  The  analysis  is 
based  on  linear  elastic  response,  and  the  reinforcing  rods  are  considered 
to  have  different  elastic  constants  and  different  density  from  those  of 
the  matrix.  Attenuation  of  the  stress  wave  then  arises  from  refraction 
and  reflection  and  the  associated  generation  of  shear  waves-  The  shape 
of  the  front  ceases  to  be  plane,  and  this,  as  well  as  the  attenuation, 
influences  the  tensile  stress  magnitudes  caused  by  reflection  at  the 
back  free  surface.  The  analysis  presented  here  applies  to  general 
three-dimensional  wave  propagation  with  inclusions  of  arbitrary  shapes, 
although  the  examples  given  involve  only  two  space  variables  for 
illustrative  purposes. 

We  consider  an  initial  wave  as  shown  in  Fig.  lb  with  a  discon¬ 
tinuous  rise  in  stress  normal  to  the  wave  front  of  magnitude  oQ 
from  the  unstressed  state  ahead  of  the  wave.  Behind  the  front 
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the  stress  is  considered  to  decrease  gradually.  Whether  subsequently 
the  stresses  at  the  wave  front  are  critical  from  the  damage  standpoint 
depends  on  the  detailed  structure  of  the  composite  medium,  the  properties 
of  its  components,  and  the  width  of  the  slab.  It  can  happen  that  redis¬ 
tribution  of  the  energy  in  the  wave  leads  to  concentration  of  stress 
well  behind  the  wave  front,  and  the  present  approach  would  then  no z 
elucidate  the  dominant  mode  of  fnil’vre.  The  present  solution  indicates 
when  the  wave  front  will  be  suffics "ntly  attenuated  for  this  circumstance. 

The  basis  for  the  present  ap1  oach  is  the  ray  tracing  and  associated 
wave  front  analysis  of  geometric  optics  [l~5]*  It  can  be  shown  that 
even  in  the  case  of  elastic  waves  with  two  basic  wave  speeds  in  each 
material,  the  stress  amplitude  at  the  initial  wave  front  is  given 
exactly  by  the  laws  of  geometrical  optics  associated  with  the  propaga¬ 
tion  of  irrotational  waves.  This  follows  from  application  of  the  theory 
of  characteristics  and  is  in  accord  with  the  more  common  interpretation 
of  geometrical  optics  as  an  approximation  to  the  solution  for  oscillatory 
waves  at  high  frequencies.  This  approach  determines  the  stress  tensor 
at  the  wave  front,  and  can  be  extended  also  to  determine  stress  gradients 
there . 

It  is  well  known  that  the  displacement  vector  u  for  dynamic 
analysis  of  an  isotropic  linear  elastic  body  can  be  expressed  in  terms 
of  the  scalar  potential  ®  and  the  vector  potential  £  in  the  form 

u  =  VG>  +  VX£  (l) 
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where  <D  and  £  satisfy  the  following  equations 


V 


1 


Q2  =  ( X+2  |i)/  p 


(2) 


P2  =  lV  P 


(3) 


V  £  =  0 

rs/ 


(4) 


X  and  p.  are  Lame's  constants  and  P  is  the  density  of  the  material. 
Since  cc  >  p,  the  scalar  potential  $  is  governed  by  a  wave  equation 
with  a  greater  wave  velocity,  and  hence  for  a  wave  moving  into  an 
undisturbed  medium,  the  region  adjacent  to  the  leading  wave  front 
corresponds  to  zero  £,  until  the  slower  distortion  wave  arrives. 
This  follows  from  the  "region  of  influence"  theory  of  the  solutions  of 
hyperbolic  equations  [6].  Therefore  u  =V<I>  behind  the  wave  front. 

Let  (1^,  be  a  system  of  orthogonal  coordinates  with 

normal  to  the  wave  fronts.  Thus  =  constant  is  the  front  of  the 
disturbance.  If  we  consider  $  to  be  zero  in  the  undisturbed  region 
ahead  of  the  front,  0  is  also  zero  immediately  behind  the  front. 

Since  displacements  must  be  continuous  in  an  elastic  body, 
u  =v$  =  0  immediately  behind  the  wave  front.  This,  together  with 
the  fact  that  <1>  =  0  at  the  wave  front  gives 


d24> 


except  for 


i=j=l  • 


1+ 


Accordingly,  the  only  non-vanishing  strain  component  immediately  behind 
the  wave  front  is  e^.  For  this  simple  one-dimensional  strain,  the 
stresses  can  be  written  down  immediately  using  Hook' s  law: 


"'ll  =  (1+2>*)£11  »  °22  =  °33  ■  Xcll 


"12  *  "23  ’  "15  '  0  ' 


Introducing  the  average  hydrostatic  tension 


(5) 


p  =  aii/3  =  (X  +  I  |1)eii  =(X+|  p)?  «  , 


(6) 


since  in  the  present  case,  Eq.  (5)  is  reduced  no. 


_  _  X+2p 

11  5x^2p  3p  * 


°22  =  °33  "  3>,+2p 


5p 


°12  =  °25  =  °33  =  ° 


(7) 


p 

Thus  it  suffices  to  study  the  propagation  of  p,  or  v  4  by  Eq.  (6), 
which  satisfies  the  acoustic  wave  equation  (2),  in  order  to  determine 
the  variation  of  stress  generated  by  the  passage  of  the  wave  front. 

The  problem  therefore  is  reduced  to  acoustic  wave  analysis  with  a 
single  wave  velocity  as  for  a  gas,  and  the  conditions  immediately  be¬ 
hind  the  wave  front  can  be  determined  by  geometrical  acoustics  for  the 
single  scalar  wave  equation. 
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2. 


Geometrical  Acoustics . 


Ray  tracing  and  wave -front  analysis  for  the  scalar  wave  equation 
has  been  termed  "geometrical  acoustics"  by  Priedlander  [2]  in  a  detailed 
presentation  of  the  topic.  He  developes  the  theory  of  characteristics 
and  shows  that  wave  fronts  as  they  vary  in  time  form  a  three-dimensional 
characteristic  manifold  in  four  dimensional  space-time.  If  t  =  t(x) 
is  the  equation  for  the  wave  front  at  time  t,  it  is  known  that  f(x) 
satisfies  the  eikonal  equation 

(VT)  •  (VT)  =  (8) 

a 

The  rays  normal  to  these  fronts  are  bi -characteristics,  and  a  system 
of  ordinary  differential  equations  governs  conditions  along  each  ray. 

The  propagation  of  discontinuities  is  considered  through  the  analysis 
of  "weak  solutions"  of  the  wave  equation  using  the  theory  of  distribu¬ 
tions  or  generalized  functions  [?]•  It  is  shown  that  the  saltus  or 
jump  across  a  discontinuity  satisfies: 

2(Vt)  •  V[p]  +  [p]^x  =  0  (9) 

which  Luneburg  [l]  called  the  transport  equation,  [p]  stands  for  the 
saltus  which  is  identical  to  p  behind  the  wave  front  in  the  present 
case,  since  p  is  zero  ahead  of  the  wave  front.  The  transport 
equation  (9)  can  be  reduced  to  an  ordinary  differential  equation  along 
the  ray  and  yields,  when  integrated,  the  intensity  law  of  geometrical 
optics,  that  the  energy  density  integrated  across  a  tube  of  rays 


6 


remains  constant.  For  propagation  through  a  homogeneous  medium,  the 
energy  can  be  represented  by  the  square  of  the  hydrostatic  tension  p, 
so  that  p  varies  inversely  as  the  square  root  of  the  area  of  cross 
section  of  the  ray  tube.  For  a  homogeneous  medium  the  rays  are  straight, 
so  that  the  area  of  an  infinitesimal  ray  tube  is  given  directly  in  terms 
of  the  Gaussian-  curvature  of  the  wave  front.  If  R  and  S  are  the 
principal  radii  of  curvature  at  the  wave  front,  we  have 

p  RS  =  constant  (10a) 

along  the  ray.  If  one  of  the  radii,  say  S,  is  infinite  as  in  the 
case  of  two-dimensional  wave  motions,  we  have,  instead  of  (10a) 

2 

p  R  =  constant  (10b) 

Changes  in  intensity  at  the  boundary  between  the  matrix  and  inclusions 
are  determined  by  the  equivalent  of  Fresnel' s  formulae  for  elastic 
media  [7] • 

The  above  rules  are  identical  with  those  used  in  classical  geome¬ 
trical  optics  and  considered  as  an  approximation  to  the  propagation 
of  high  frequency  waves.  It  was,  we  believe,  first  pointed  out  by 
Luneburg  in  the  lecture  notes  which  comprised  the  fore-runner  of  [l], 
that  they  express  the  propagation  of  a  discontinuity  exactly  within 
the  limits  of  linear  theory.  In  che  following  two  sections,  we  illus¬ 
trate  the  method  of  geometrical  acoustics  by  considering  a  plane  sharp 
wave  front  of  the  type  shown  in  Fig.  lb  propagating  through  an  elastic 
circular  cylindrical  inclusion.  The  case  when  the  inclusion  is  a  sphere 
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is  discussed  in  Section  5 •  In  Section  6,  we  derive  the  equations  which 
enable  one  to  determine  the  stress  gradient  and  the  higher  order  deriva¬ 
tives  of  stresses  behind  the  wave  front. 

3 •  Geometry  of  the  Wave  Fronts. 

Fig.  2  shows  the  wave  front  configurations  generated  when  the  plane 
incident  wave  interacts  with  an  embedded  reinforcing  circular  cylindrical 
inclusion.  The  ratio  Q^/o^  of  the  dilatation  wave  speed  in  the  rein¬ 
forcement  (medium  2)  to  that  in  the  matrix  (medium  l)  is  taken  to  be 
4/3*  As  the  wave  front  strikes  the  inclusion  both  refracted  and  reflected 
<E>  and  £  waves  are  produced.  As  shown  in  Fig.  2a  an  incident  wave  thus 
produces  four  different  wave  fronts.  The  motion  in  the  regions  between 
the  dilatation  wave  front  (<J>)  and  the  distortion  wave  front  (£) 
is  irrotational  because  £  is  zero.  However,  this  does  not  imply 
*hat  shear  stresses  are  zero-  In  view  of  symmetry,  only  upper-halves 
of  the  wave  fronts  are  shown  in  the  figure. 

The  wave  fronts  become  more  complicated  when  the  incident  wave 
front  passes  the  point  Q  where  total  reflection  of  the  dilatation 
wave  occurs  and  the  incident  wave  produces  no  refracted  dilatation  wave 
(Fig.  2b).  For  simplicity  distortion  wave  fronts  are  omitted  from  this 
figure.  Also  omitted  are  the  wave  fronts  due  to  multiple  refractions 
and  reflections  except  the  precursor  dilatation  wave  abed  in  Figs. 

2b  and  2c.  All  rays  which  reach  the  wave  front  cd  in  Fig.  2b  have 
undergone  a  sing.le  refraction,  while  the  rays  which  reach  the  wave  front 
ebc  have  been  subjected  to  two  refractions  through  the  inclusion. 
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In  the  following  our  main  interest  will  be  the  determination  of 
the  wave  front  abed  and  the  discontinuity  in  pressure  at  this  wave 
front.  Therefore,  we  need  consider  only  those  incident  rays  which 
meet  the  inclusion  at  points  below  Q.  The  reflected  wave  front  ef 
can  be  deteirainea  in  the  same  manner. 

Although  the  wave  fronts  can  be  located  by  integrating  Eq.  (8),  it 
is  simpler  by  a  direct  calculation  using  Fermat's  principle.  To  illus¬ 
trate  the  procedure  in  computing  rays  and  wave  fronts,  let  us  consider 
a  ray  AB  in  Fig.  3  which  strikes  the  inclusion  at  angular  position 
0.  If  a  is  the  radius  of  the  inclusion,  the  coordinates  of  B  are 
given  by: 


Xg  =  a(l-cos9) 

(11a) 

yB  =  a  sin® 

(lib) 

Let  t  =  0  be  the  time  at  which  the  initial  plane  wave  reaches  x  =  0 
and  contacts  the  inclusion.  Then  the  time  t  required  for  point  C 
on  the  refracted  ray  BD  to  be  reached  by  the  wave  front  is 


t  = 


+bo 

°i  “a 


or 


V  -  t  -  S >  (12> 

where  ^  is  the  length  of  the  straight  ray  BC.  The  inclination  of 
the  refracted  ray  to  the  incident  ray  direction,  5,  is  given  by  Snell' s 
law: 
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(15) 


sin(9  +  5)  _  ^2 
sin9  =  0^ 

The  coordinates  of  the  point  C  are  thus: 

XC  =  ^  +  Vos  6  (1^a) 

=  a(l-cos0)  +a  cos&[t  -  (l-cos0)3 

d  “l 

yC  =  yB  +  £BCSin5 

=  a  sin©  +  a  sin  8  [t  -  —  (1  -  cos©)] 

1 

For  a  given  t,  (l4a,b)  is  the  equation  for  the  wave  front  cd  in 
Fig.  2b  with  0  as  parameter.  By  direct  calculation  it  can  be  shown 
that  the  slope  of  the  front  is  perpendicular  to  the  ray  BC,  as  it 
must  be  according  to  Huygens'  construction  for  wave  fronts. 

A  similar  calculation  determines  the  wave  front  after  refraction 
out  of  the  inclusion.  For  the  coordinates  of  point  E  in  Fig.  j>,  we 
have 

XE  =  XB  +  i!BDCOS&+^DECOs25  (l5a) 

yE  =  yB+'eBDSin6+iDESin26  (l5b) 

where 

ilBD  =  2a  cos(0+8)  (l6) 

and  i__  is  determined  by 
lib 
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or 


0!, 


*DE  =  °lt 


XB  "  a2  BD 


(17) 


Clearly,  one  can  continue  this  procedure  to  determine  the  wave 
fronts  if  there  is  more  than  one  inclusion. 

The  wave  fronts  as  expressed  by  Eqs.  (l4a,b)  and  (l5a,b)  are  shown 
in  Fig.  3  for  various  times.  The  dotted  line  Qb  is  the  locus  of 
intersection  of  the  initial  incident  plane  wave  front  and  the  refracted 
wave  front.  The  determination  of  the  pressure  p  on  the  wave  front 
will  be  described  in  the  next  section. 


4.  Stress  Magnitudes  at  the  Wave  Front. 

The  change  in  magnitude  of  stresses  due  to  refraction  at  point 
B  or  D  in  Fig.  3  should  be  independent  of  the  curvatures  of  the 
incoming  wave  front  and  the  interface.  While  this  statement  seems 
intuitively  correct,  there  appear  to  be  no  proofs  available.  We  will 
prove  this  in  Section  6,  and  therefore  we  accept  this  fact  as  valid  in 
this  Section. 

For  a  given  plane  dilatation  wave  ^  of  magnitude  A.^  incident 

on  a  plane  interface,  the  refracted  dilatation  wave  O'  of  magnitude 

A'  has  been  determined  in  [7].  However,  we  are  interested  in  p,  and 
2 

hence  v  $  by  Eq.  (6).  If 
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£L  £L 

*  "A1 


at  the  interface,  then  according  to  [73, 


2^.  =  410?  (18) 

^  A^a' 

Making  use  of  Eq.  (6)  and  the  result  obtained  in  [73  for  the  value  A' /A^, 
we  write  the  pressure  ratio  p’ / p  at  the  interface  in  the  following 
form: 


f  -  p<e’ 


a 


&  £ 
a  ;  a' 


p'(a:'a  -ip'V  /2cos9M\ 

"  P(a2  -|rV2  I™*1*'1  3 


(19) 


K  =  cos  9  (~  +  r^sin  S)  +  ^r 


Jl-(~  sin 9)2  (l-qsin2e) 


M  =  (~  +  qsin20)  + 


I i-(  j-  sine  )2 
J 1-(|  sine)2 


(l-qsin20) • 


L  =  §  - -  (i  .  £l-^sin20)  -  |j-  qsine  fl-(~  sine)2 

Jl-(g  sine)2 

N  =  (l  -  — ■ - qsin2e)sine  -  rjsinecose  sin e)2 


where 
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4f>2  -  f 


The  quantities  with  a  prime  refer  to  the  refracted  medium  while  the 
quantities  without  a  prime  refer  to  the  incident  medium.  The  trans¬ 
mission  coefficient  as  expressed  by  Eq.  (19)  has  been  tabulated  in 
[8]  for  —  =  =  1//J  with  variable  Q,  and  . 

Let  (Pg)2  and  (Pg)^  be  pressure  at  point  B  of  Pig.  3 
on  the  side  of  medium  2  and  medium  1  respectively.  Since 

( p_ )  =  p  ,  the  incident  constant  pressure,  (p  )  can  then  be 

B  ±  O  B  d 

determined  by  using  Eq.  (19).  The  geometrical  acoustics  law  Eq.  (10b) 
then  gives 


where  (Kg)2  is  the  radius  of  curvature  of  the  wave  front  at  the 
point  B  in  medium  2,  and  R  is  the  radius  of  curvature  at  0. 

V-» 

Because  the  rays  are  straight  lines  normal  to  the  wave  front  in  each 
medium 


EC  ■  <V2 +  V  (21) 

and  hence  we  can  determine  the  pressure  at  any  point  in  the  inclusion. 
Although  (Rp)p  could  be  obtained  by  differsntiation  from  the  wave 
front  expression  (l4),  it  is  more  convenient  to  analyse  the  refraction 
geometry  adjacent  to  the  point  B. 

Let  us  consider  a  more  general  situation  as  shown  in  Fig.  4  in 
which  an  incident  ray  OB  strikes  the  interface  boundary  at  B  and 
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produces  a  refracted  ray  BC.  Let  fL  and  denote  respectively  the 
radii  of  curvature  of  the  incident  and  refracted  wave  fronts  at  point 
B.  Moreover  let  denote  the  radius  of  curvature  of  the  interface 

boundary  at  point  B.  By  considering  an  adjacent  incident  ray  OD, 
it  can  be  shown  that  the  following  relation  holds: 

COS0 


cot  0. 


/COS0.  , 
i+1 


R. 


"b 


=  cot  9 


R 


r  +  1 


*b 


(22) 


In  Eq.  (22),  R^  assumes  a  negative  value  if  the  center  of  curvature 
of  the  interface  is  on  the  left  of  the  boundary  instead  of  on  the  right 
as  shown  in  Fig.  4.  Likewise,  R^  and  R^  can  assume  negative  values. 

Applying  Eq.  (22)  to  point  B  in  Fig.  3,  we  obtain,  since 
R^  =  co,  and  using  Snell's  law,  Eq.  (13), 


<Va  -  a*  eos2(8*6>  (25) 

cos0-cos(0+6) 

1 

With  Eqs .  (2l)  and  (23),  Eq.  (20)  determines  the  pressure  p  at  the 

wave  front  inside  the  inclusion.  Figure  5  shows  an  example  of  the 

pressure  p  as  a  function  of  position  when  the  wave  front  arrives* 

The  results  are  shown  in  the  form  of  contour  lines  of  constant  pressure. 

The  material  constants  used  are  a^/a^  =  ~  $2/^2.  ~  and. 

Pg/p^  =  1.  The  pressure  at  point  0  inside  the  inclusion  is  1.142  pQ 

and  increases  to  2.522  pQ  as  Q  is  approached  from  0  along  the 

circumference  of  the  inclusion.  On  the  other  hand,  the  pressure  at 

H  inside  the  inclusion  is  O.885  p  and  decreases  to  zero  as  we 

o 

approach  Q  along  the  circumference.  Thus  Q  is  a  point  of  singularity. 
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Figure  5  also  shows  the  results  of  similar  calculations  for  pressure  at 
the  wave  front  after  refraction  out  of  the  inclusion.  The  pressure  at 
H  outside  of  the  inclusion  is  0.758  pQ  and  decreases  to  zero  as  we 
approach  Q  along  the  outside  of  the  circumference.  The  line  Qe  is 
the  contour  line  for  zero  pressure.  However,  if  we  regard  Qe  as  the 
reflected  ray  at  Q  and  calculate  the  reflected  pressure  along  Qe,  we 
will  find  that  the  reflected  pressure  along  Qe  is  not  zero.  Hence 
there  is  a  discontinuity  in  pressure  across  Qe.  In  fact,  the  radius 

of  curvature  of  the  wave  front  feb  of  Fig.  2b  which  consists  of  the 

reflected  front  fe  and  the  refracted  front  eb  is  discontinuous  at  e. 

Figure  5>  when  used  in  conjunction  with  Fig.  3>  enables  one  to 

determine  the  pressure  at  the  wave  front  for  a  given  time  t.  The  wave 
fronts  obtained  in  Fig.  3  are  also  shown  in  Fig.  5.  Notice  that  the 
pressure  in  the  region  above  the  dotted  line  Qb  is  the  pressure  on 
the  refracted  wave  front.  The  actual  pressure  should  include  the 
incident  plane  wave  pressure  pQ.  The  stresses  at  the  wave  front  in 
this  region  should  be  modified  accordingly. 

The  pressure  along  the  x-axis  has  a  simple  form  and  is  recorded 
here.  Substituting  0=5=0  in  Eq.  (23),  we  have 


R  + 
0^ 


(2k) 


where  the  o+  refers  to  the  abscissa  value  of  the  point  designated. 
Evaluation  of  (19)  gives 


-  .  ■ - - - 
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V 


1 


1 


jJfi 


and  application  of  (20)  yields 


(25) 


p(x)  =  (p  x) 


a  +  x(~  -  1) 

1 


,  0  <  x  <  2a 


Refraction  out  of  the  inclusion  gives 


R2a-  =  2R2a+ 


(26) 


(27) 


and  hence 


p(x)  = 


U  +  f£i  ♦  /(a^-D  tatS-aj^-D 

\  »ift>  »iaJ'/  “i  *  “i 


,  x  >  2a 


(28) 


until  the  next  reinforcing  cylinder  is  reached.  It  is  interesting  to 
note  that  for  O.g  >  Ct^,  the  usual  case  for  composites,  the  wave  fronts 
diverge  and  the  stress  amplitude  attenuates.  However,  for  ag  <  0^,  a 
converging  wave  front  occurs  after  the  first  refraction,  and  the  energy 
becomes  concentrated  at  the  caustic  giving  infinite  stress  values  accord¬ 
ing  to  the  linear  theory.  Stress  concentration  can  occur  also  at  the 
distortion  wave  front  when  fig  <  Ct  • 


5 •  Wave  Propagation  through  a  Spherical  Inclusion. 

For  an  initial  plane  wave  as  shown  in  Fig.  lb  propagating  through 
a  spherical  inclusion  of  radius  a,  the  geometry  of  the  wave  fronts 
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intersecting  a  plane  through  the  axis  of  symmetry  is  the  same  as  for  the 
cylindrical  inclusion  shown  in  Pig.  3*  Thus  the  results  obtained  in 
Section  3  apply  equally  well  to  a  spherical  inclusion  if  we  replace  y 
by  r,  the  radial  distance  from  the  axis  of  symmetry  x.  The  proce¬ 
dures  outlined  in  Section  4  for  determining  the  stress  magnitude  at 
the  wave  front  are  still  applicable  except  that  the  intensity  law 
Eq.  (10a)  should  be  used  instead  of  Eq.  (10b).  Thus,  Eq.  (20)  should 
be  replaced  by 


PC  _  B W2" 
(pB^2  v  RCSC 


(29) 


where 


(sb)2 


+  £. 


BC 


(30) 


(a  )  -  a  sing 

VV2  “  sin 6 


(31) 


In  Fig.  6,  we  show  the  pressure  at  the  wave  front  for  a  spherjcal  inclu¬ 
sion  using  the  same  material  parameters  as  in  Fig.  5  for  comparison. 
Again,  point  Q  is  a  point  of  singularity  and  the  pressure  suffers  a 
jump  across  Qe.  As  is  to  be  expected,  more  rapid  attenuation  of  the 

wave  front  occurs  for  a  spherical  inclusion  than  for  a  cylindrical 
inclusion  because  of  spreading  in  all  directions  parallel  to  the  wave 
front,  and  the  formulae  for  attenuation  along  the  x-axis  Eq.  (26,28) 
are  simply  modified  by  removal  of  the  square  root  since  to  second 
order  the  wave  front  is  spherical  by  symmetry  on  this  axis* 
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2£ 


mim 


a-tr 


6.  Stress  Gradient  and.  Higher  Order  Derivatives  Behind  the 

Wave  Front. 

To  determine  the  spatial  derivatives  of  stress  right  behind  the 
wave  front,  we  have  to  establish  the  higher  order  transport  equations 
which  determine  the  propagation  of  discontinuities  in  stress  derivatives 
along  the  ray  and  an  extension  of  the  equivalent  of  Fresnel' s  formulae 
for  the  reflection  and  transmission  coefficients  of  the  stress  deriva¬ 
tives  at  an  interface.  It  should  be  noted  that  spacial  derivatives  and 
time  derivatives  of  any  quantity,  say  3>,  behind  the  wave  front  are 
related  once  we  know  ®  at  the  wave  front.  For,  if  t  =  i(x)  is  the 
wave  front  at  time  t,  we  have 

V($>(x,t(x)}  =  tv<I>(x,t)}t=T(x)  +Vt(x)  *(2S,t))t=T(s)  (32) 

Similarly,  the  k-th  order  spatial  derivative  of  G>  at  the  wave  front, 
the  limits  being  taken  as  the  front  is  approached  from  behind  which 
applies  throughout  this  section,  can  be  determined  once  we  know  the 
first  k  order  time  derivatives  of  $  at  the  wave  front.  It  suffices 
therefore  to  discuss  the  propagation  of  discontinuities  in  time  deriva¬ 
tives  • 

Let  t  =  t(jj)  and  t  =  t(x)  be  the  dilatation  and  distortion 
wave  fronts  respectively  at  time  t.  Then  t(x)  and  ”(x)  satisfy 
the  eikonal  equations 

(vt)  *  (Vt)  =  ~  (35a) 

Q 

(vT)  •  (?7)  =  -L  (33*>) 

r 
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Define 


i  =  t  -  ?(x) 

(3*0 

r-t  -  t(x) 

(54b) 

Physically, 

t ’  and  £  are  the  times  elapsed  after  the  arrival  of  the 

dilation  and 

distortion  wave  fronts.  At  the  dilatation  wave  front 

i  =  0  and  at  the  distortion  wave  front  £  =  0.  Now  we  expand  4> 

and 

J  into  the  following  series  [2]: 

rw 

<K&t)  =  £  q>k(x)  i’k+2(0 

(55a) 

jk(»)  fkta(D 

(55b) 

where  fQ(o)) 

is  an  arbitrary  function  of  cd  for  a>  >  0  and  vanishes 

for  w  <  0, 

and 

r(ti 

fk+1(m)  =  /  fk(a>')dco’ 

J  -00 

(56a) 

or 

—  f.  (o)  =  f,  (a>) 
dn>  k+1'  '  k  ' 

(56b) 

Thus  fk(o)  also  vanishes  for  o  <  0  for  all  k  >  0.  Notice  also 
that  fk(o)  =  0  for  all  k  >  1.  In  Eqs.  (35)  the  terras  involving  f 
and  f^  are  excluded  because  of  the  continuity  requirement  for  dis¬ 
placements  across  a  wave  front. 
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The  expansion  expressed  by  Eqs*  (35)  has  a  physical  meaning  if  we 
choose  f  (co)  =  H(co),  the  Heaviside  step  function.  For  then  f,  (us)  = 

O  K 

a/  ki  and  q>k(j$)  becomes  the  k-th  time  derivative  of  «1>( jc, t )  at  tne 
>ave  front  t  =  r(x).  Similarly,  i£k(x)  becomeb  the  k-th  time  deriva¬ 
tive  of  i(x,t)  at  the  wave  front  t  =  t(x). 

By  substituting  Eq,  (35a)  into  Eq,  (2)  and  making  use  of  Eqs. 
(33a)  and  (36b),  we  obtain  the  higher  order  transport  equation  for  cpk 
in  a  recurrence  form: 

2(Vt)  •  (^<pk)  +  q>k  72t  =  V2(Pk_1  >  (k=0, 1,2, ... )  (37a) 

provided  we  define  <p  ^  =  0.  When  k  =  0,  Eq.  (37a)  reduces  to  the 
zero  order  transport  equation  (9)-  Similarly,  we  obtain  the  higher 
order  transport  equations  for  ^k  which  can  be  written  in  a  dyadic 
notation  as 

2(V7)  *  (v£k)  +  £k  V2t  =  ,  (k=0, 1,2,  . . . )  (37b) 

provided  we  define  ,  =  0  . 

To  obtain  the  reflection  and  transmission  coefficients  at  an 
interface,  ve  have  to  find  the  displacements  and  stresses  at  the  inter¬ 
face  and  apply  the  continuity  condition.  Substituting  Eqs.  (35)  into 
Eq.  (l)  we  obtain 

b  ■  J0  {(vvr4  VT)iW£  > +  <v  *  4-f v7x  4>  WD)  (38) 

Since  t,  =  £  at  the  interface,  we  have 
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(39) 


Ub  -  i  {-(fkvT  +  vt 

k=0 


where 


4-1  “  VVl  +7X  4-1 

which  contains  only  the  lower  order  terms  and  a  superscript  b  stands 
for  the  interface  boundary.  If  we  substitute  Eq.  (j58)  into  the  following 
stress-displacement  relation 


ff  =  )i(vu  +  u v)  +11  V  ‘  u 


(4o) 


and  let  £  =  we  oboain  stresses  at  the  boundary 


h  r  >  —  —  __  _ 

2  =  I  (2^  VT  VT  +  ~  l)(p.  +  |l(VT  VT  X  i  +  VT  X  t  Vt) 

k=o  ad  "  k  ~K  ~K 


§k-i}4«> 


(41) 


provided  we  define  <p_,,  =  0,  =  0.  contains  only  the  lower 

order  terms  ^k-2'  ~k-l'  anc^  ^k-2’  ~  an<*  ~  a 

vector  normal  to  the  dilataticn  wave  front  and  the  jdistortion  wave  front 
respectively.  Then,  by  Eqs.  (35); 


"  §  & 


(42a) 


(42b) 


Figure  7  shows  a  typical  configuration  in  which  reflected  and  transmitted 
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<5 *i1)B(1)  -  f  H(i)  x  *ii)) +  (5  4r)o(r)  +  p  ®(r)  x  «:lr)> 


<sMt)a(t)  +Fa(t)x^t))  -!S-1 


u*  =  u^  j  -  uj.1]  -  u[r] 

2.  ~k-l  ~k-l 


The  condition  that  the  surface  traction  crb*  continuous  at 

s 

point  B  gives,  using  Eqs.  (4l)  ~  (44)  and  Fig.  7, 


—  (2pcos0n(i)+Xnb)<p[i)+p(-o[i)sin0£(;L)+cosen(i)  x 

+  ~  ( -2pcos  0  n^r^+Xnb)<p^r^+p(-Q^r^sin  0n^r^-cos0  n  ^  x  V*^) 
=  -^2  (2p*  cos0'  n^b ^+X'  nb J^j^+P1  ( -oj^sin0'  n^b ^+cos0! £^b^X  ^j^^k-l 


where 


%  -  &  • 


m*  _  m(4)  _  rp(  i )  _  m(^') 

2_  rv'&  "  1  -  2. 


T,  =  S,  .  .  n 
~k-l  =jk-1  ~ 


If  we  substitute  Eq.  (35b)  into  Eq.  (4)  and  make  use  of  Eq.  (36b),  we 
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(49a) 


P  V  * 

fv 


(r) 

k-1 


6'  V 


Lk-1 


(49b) 


(r) 

Equations  (46)  ~  (49)  give  a  recurrence  formulae  for  determining  , 

iw 

<pj^,  and  (k=0, 1,2,  . .  • )  at  the  interface.  Thus  one  can 

calculate,  in  principle,  the  reflection  and  transmission  coefficients 
for  any  order  of  time  derivatives  of  stress  behind  the  wave  front. 
Equations  (46)  and  (47)  can  be  reduced  to  scalar  equations.  If 

Id 

we  perform  the  scalar  product  on  both  sides  of  Eq.  (46)  by  t  and  n 
respectively,  we  have 


((pki)+<pkr))  +  c°t0(4i)-o[r))  =  <p{Ho3te’ 


_i* _  ij*  •  fe 

sin0  ~k-l  ~ 


(50) 


eot0(9^X  )_(Qk  ^ 


)  =  ccte’4t)-0<t) 


[3  * 

— —  u.  . 

sin0 


b 

n 


(51) 


Similarly,  if  we  perform  the  scalar  product  on  both  sides  of  Eq.  (47) 
by  nb  and  t  respectively,  we  obtain: 

p{-(cot20-l)(9[i)^r))+2  cot^Q^-Q^)} 

=  p'|-(cot20'  -1)9^^+  2cot  0’Q^}  -  — ^2=  T*_x  •  nb  (52) 

sin  0 

p{2cot  0(9^  -?[r) )+( cot20 -1)  ( )+Q* r } ) } 

=  p'  (2cot  0'  9^  +  (cot2©'  -1)0^}  +  -^2=  T*  •  t  (55) 

sin  0 
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Finally,  we  perform  the  scalar  product  on  both  sides  of  Eqs,  (46)  and 
(47)  by  b  to  obtain: 


P(i) 

k 

P(i) 

k 


+  4r)  -  r  e  4-i '  6 


-  p<r)  =  eLssdL  P(t)  .  u*  ,b 

K  p  cos  0  pcos  0 


(54) 

(55) 


where 


Pk  =  (n  x  b)  •  £k  (56) 

Equations  (49)  ~  (55)  are  alternate  equations  for  (pj^,  <Pk^, 

(t) 

and  ^k  '  -  Notice  that  Qk  and  ?k  are  the  components  of  ^  in 
the  b  and  n  x  b  directions  respectively  while  Eqs*  (49)  give  the 
component  of  £k  in  the  n  direction*  Thus  we  have  three  components 
of  in  three  mutually  orthogonal  directions. 

It  is  known  that  one  can  delete  any  one  of  the  three  components 
of  without  loss  of  generality  by  ignoring  the  condition  (4).  In 
that  case,  Eqs*  (49)  which  are  derived  from  Eq.  (4)  are  not  needed  and 
Eqs*  (50)  ~  (55)  alone  are  sufficient  for  determining  the  unknowns. 

For  two-dimensional  wave  motion  and  axially  symmetric  wave  motion, 
only  one  component  of 
Equations  (54)  and  (55)  are  automatically  satisfied  and  Eqs.  (50)  ~  (55) 
alone  are  sufficient  for  determining  the  unknowns. 

It  should  be  mentioned  that  Eqs.  (50)  ~  (53)  are  identical  to 
Eqs*  (3-11)  ~  (3-14)  of  (7]  for  plane  waves  incident  on  plane  boundary, 
except  for  the  lower  order  terms  which  are  zero  when  k  =  0  and  an 


,  namely  as  defined  by  (48),  is  needed. 
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obvious  difference  in  the  use  of  notations.  Thus  the  assertion  that 
the  reflection  and  transmission  coefficients  for  the  stress  amplitudes 
at  an  interface  are  independent  of  the  curvatures  of  the  incident  wave 
front  and  the  interface  surface  has  been  proved.  The  reflection  and 
transmission  coefficients  for  the  stress  gradient  and  the  higher  order 
derivatives  in  general  depend  on  the  curvatures  of  the  incident  wave 
front  and  the  interface  surface  through  the  lower  order  terms  in  Eqs • 

(49)  ~  (55). 

7 •  Discussion. 

The  analysis  presented  here  applies  to  inclusions  of  general 
geometry  and  can  be  used  to  determine  the  stress  behind  a  wave  front 
which  is  reflected  from  an  interface.,  a  rigid  surface,  or  a  free  surface. 
For  the  latter,  the  stress  due  to  the  incident  wave  should  be  superposed 
to  obtain  the  actual  stress  behind  the  reflected  wave  front. 

The  analysis  can  also  be  applied  to  determine  the  discontinuity  in 
shear  stress  across  a  distortion  wave  front.  However,  the  geometry  of 
the  distortion  wave  front  is  in  general  more  complicated  than  the  dila¬ 
tation  wave  front.  Figure  8  shows  some  of  the  distortion  wave  fronts 
for  the  same  material  parameters  as  in  Figs.  3,  5 >  and  6.  For  compari¬ 
son,  the  dilatation  wave  fronts  and  the  corresponding  rays  obtained  in 
Fig.  3  are  shown  as  dotted  lines.  P  is  the  point  at  which  the  incident 
ray  is  tangential  to  the  circumference  of  the  inclusion.  After  passing 
point  Q,  there  are  two  distortion  wave  fronts,  one  is  formed  by  one 
refraction  at  the  interface  while  the  other  is  formed  by  reflection  of 
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the  refracted  dilatation  wave  front  inside  the  inclusion.  The  curve 
qs  is  the  caustic  of  the  reflected  distortion  rays  and  therefore  shear 
stress  becomes  infinite  along  qs •  The  shear  stress  is  also  infinite 
at  the  wave  front  where  the  associated  rays  have  touched  the  caustic 
qs.  (see  [2,3]).  This  situation  does  not  seem  to  have  been  mentioned 
previously  in  the  literature  on  composite  materials.  In  fact  if  one 
considers  the  second  dilatation  wave  front  which  is  obtained  from  the 
reflection  of  the  refracted  wave  front  inside  the  inclusion,  one  again 
obtains  a  caustic  similar  to  the  curve  qs  in  Fig.  8.  Thus  we  also 
have  an  infinite  stress  in  the  inclusion  due  to  dilatation  waves.  For 
a  real  material,  this  implies  either  a  failure  in  the  inclusion  or 
the  material  near  the  curve  qs  becomes  plastic  and  hence  subsequent 
wave  propagation  analysis  should  be  modified.  In  experiments,  an 
infinite  stress  concentration  along  the  curve  qs  is  unlikely  to 
happen  because  it  is  not  possible  to  generate  a  truely  sharp  wave  front 
as  shown  in  Fig.  1(b).  If  there  is  a  finite  time  for  the  stress  to 
rise  to  its  peak  magnitude,  an  infinite  stress  will  not  occur  along  qs. 
However,  this  does  not  rule  out  the  possibility  of  an  infinite  stress 
concentration  occurring  along  other  caustics  which  are  formed  by  multiple 
reflections  inside  the  inclusion.  For  instance,  an  incident  pulse  with 
initially  continuous  stress  at  the  wave  front  can  become,  after  the 
first  reflection  inside  the  inclusion,  a  shock  wave  front  with  a  finite 
discontinuity  in  stress  at  the  wave  front.  This  shock  wave  front  can 
then  become  a  front  with  an  infinite  stress  after  the  second  reflection 
inside  the  inclusion. 
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It  should  be  noted  that  whether  the  stress  inside  the  inclusion 
becomes  infinite  or  not,  the  solution  presented  in  this  paper  concerning 
the  leading  wave  front  is  still  valid.  This  is  obvious  because  of  the 
hyperbolic  nature  of  the  wave  equation  considered  here.  The  discussion 
presented  above  simply  points  out  the  fact  that  a  complete  stress 
analysis  inside  an  inclusion  due  to  the  passage  of  a  pulse  is  a  com¬ 
plicated  pr  >lem.  When  an  infinite  stress  occurs,  a  realistic  solution 
can  be  obtained  only  if  one  considers  non-linear  elastic  deformation 
or  plastic  deformation. 
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WAVE  FRONT 


FIGURE  5  PRESSURE  AT  THE  WAVE  FRONT  FOR  A  CYLINDRICAL  INCLUSION 


FIGURE  7  REFLECTION  AND  REFRACTION 
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